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We calculate the current-pressure relation for pinholes connecting two volumes of bulk superfluid 
''He-B. The theory of multiple Andreev reflections, adapted from superconducting weak links, leads 
to a nonlinear dependence of the dc current on pressure bias. In arrays of pinholes one has to take 
into account oscillations of the texture at the Josephson frequency. The associated radiation of spin 
waves from the junction leads to an additional dissipative current at small biases, in agreement with 
measurements. 
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Weak links between superconducting metals have been 
extensively studied during the past 40 years. More 
recently, similar experiments have also been made in 
atomic superfluids, in the fermion '^He liquid in particu- 
lar Q . Experiments in superfluids can shed new light on 
general weak link phenomena, due to the availability of 
parameter ranges which are not easily realized in super- 
conductors. There is also intrinsic interest in superfluids 
since completely new phenomena may appear which are 
not possible in conventional superconductors. Previous 
work on ^He weak links has mostly concentrated on de- 
termining the equilibrium properties, like the Josephson 
current-phase relation I{4>). A new phenomenon iden- 
tified in this context is the anisotextural eff'ect, where 
spin-orbit degrees of freedom change as a function of the 
phase difference (j). As a consequence I{(j>) deviates essen- 
tially from the sinusoidal form I = Ic sin (jj so that there 
appear tt states, where both derivatives /'(O) and 
are positive [3, y| ■ 

In this letter we present a theory of dissipative currents 
in superfluid '^He weak links. Experimental studies of 
such time-averaged, or dc currents under a pressure bias 
have been reported in Refs.|3 and0. At constant chemi- 
cal potential difference U the supercurrent I{(f>) oscillates 
at the Josephson frequency ujj — 2U/h. Therefore the 
superfluid part does not contribute to the dc current un- 
less there is a mechanism for absorbing the extra energy 
2U per transported Cooper pair. One example is the pro- 
cess of multiple Andreev reflections (MAR), used for ex- 
plaining properties of superconducting weak links Q ■ 
Due to Andreev reflection there are quasiparticle states 
bound to the contact region. These states carry the cur- 
rent, one Cooper pair per one back-and-forth reflection. 
Simultaneously the quasiparticles gain energy until they 
escape above the gap or are relaxed by scattering. This 
process works also in superfluid '^He, but with important 
differences in details. Below we argue that in ^He there 
is an additional dissipation mechanism. Owing to the 
anisotextural effect, energy can be carried away from the 
junction by spin-wave radiation. We show that the two 
mechanisms together can explain the essential features of 
the experimental results in Refs.0 and0. 

Multiple Andreev reflections. — Let us consider a weak 



link under constant difference U in chemical potentials. 
In ^He U is related to the pressure difference P by U = 
(m^/p)P, where p is the density and ma the mass of 
an atom. The time derivative of the phase difference 
(j) across the link is given by the Josephson frequency 
(j) = LOj ~ 2U /h. The time-dependent mass current can 
then be expressed as a Fourier series 

oo 

I{P, t) = /o(P) + Y.^ll{P) sm{nujt) -f /C(P) cos(nL^jO] 

n=l 

(1) 

with coefficients loilfn In - We have calculated these co- 
efficients for a pinhole, where the dimensions of the hole 
are assumed small in comparison with the superfluid co- 
herence length ^0 = 77 nm. (All parameters are evalu- 
ated at vapor pressure.) Our calculation uses nonequilib- 
rium Green's functions, and is explained in the appendix. 
However, a large part of the results can be understood by 
the following adiabatic model, which is a generalization 
of the one used for superconductors 0- 

The current arises from quasiparticle states that are 
bound in the weak link at energies below the gap A. The 
mass current is given by 

Here eijnai'f) is the bound state energy that depends on 
the channel index i, momentum sign 5 = ±, and spin 
index cr = ±. The index n is for the case that there exist 
more than one bound state. The distribution function 
PiSnait,4') obeys the kinetic equation 

77: 1-*? a, = f /(eifag-) -p»fe<T , (3) 

at ocp 

where /(e) = [exp(e/A:BT) + l]^^ is the Fermi distribution 
and a relaxation time. 

Application of the above to ^He includes a few speciflc 
features. First, the order parameter is always suppressed 
near walls. We calculate it self-consistently, assuming 
diffusive scattering of quasiparticles at the wall. This im- 
plies that also the energies eiSna {4>) have to be calculated 
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FIG. 1: The bound state energies eiSna as a function of 
the phase difference (f). Neglecting gap suppression es = 
—5Acos{(j)/2) (dotted line). Parameters are: quasiparticle di- 
rection cosine 0.93, temperature O.QTc and R''"^ = R(=pz,&_t). 



numerically. One example is shown in Fig. ^ Second, 
since we consider the B phase, the bulk order parameter 
is of the form ARexp(i(/)), where R is a rotation matrix 
01 . Depending on the matrices on the left (1) and right 
(r) hand side&\ R*'*^, there is spin-splitting: eiSn-i4') ~ 
^iSn+i4' — "0) US' Third, the only important source 
of inelastic processes is quasiparticle-quasiparticle colli- 
sions. Near a surface the order parameter is suppressed 
and therefore we approximate the scattering rate by the 
normal Fermi-liquid form F — a[(7rfcBT)^ -I- e'^] / {tqtt^ k^) . 
Here tq — 1.14 /ismK is obtained from viscosity mea- 
surements and a is a coefficient on the order of unity. 
The low-energy F is denoted by Fq — uT'^/to. 

We are now equipped to analyze Fig.|21 Similar to the 
experiments, we study the low-bias region, U <C ksTc. In 
this region we can neglect the normal current from ener- 
gies above A, excluding only a narrow temperature slice 
near the superfluid transition temperature Tc- A charac- 
teristic scale for U is set by the scattering strength HT. 
At U = 0, the distribution function is in equilibrium, 
PiSna- = fi^iSna)- Equation ^ then gives the equilib- 
rium current-phase relation /((/i), dominated by sm((f>) 
(at least for R' = R"^) but with smaller admixtures of 
sm{n(j)) where n > 1. At U <ti hT the kinetic equation 
Q can easily be solved by linearization. This leads to a 
time-independent component Iq that is linear in U and 
independent of R^'"^, 

Io=9{T){A/hro)GnU. (4) 

Here Gn is the normal-state conductance and g{T) a fac- 
tor on the order of unity. Neglecting gap suppression, 
g{T) = J^^ tauh{Ax/2kBT){x/Vl - x^)dx. The validity 
region of the linear dependence vanishes as T ^ and 
is replaced by Iq oc C/^/^. Finally, at U HT (but still 
U <^ A) the scattering can be neglected and p is deter- 
mined by thermalization at gap edges, PiSna- — fi~SA). 
Here Iq approaches a constant value that is on the same 
order as the critical current. For a self-consistent order 
parameter all higher harmonics are effectively damped 
by the smooth eiSna-{4')- 




FIG. 2: Current amplitudes Jo, if, Jf Bind if (0 as functions 
of the chemical potential difference U. The average current 
Jo is shown at different temperatures, other curves are for 
T/Tc = 0.5. Jo is linear at small U @ and saturates at 
larger U, as discussed in the main text. The most accurate 
results are shown by solid lines. The effect of neglecting the 
gap suppression is shown by dashed lines, and the effect of 
neglecting the energy dependence of F is shown by a dotted 
line. The largest shown bias is P « 4.2 Pa. Parameters are 
a = 1.6 and R' = R". 



The Io{P) in Fig. [21 can be compared with the mea- 
surements in Ref. |j. They have the same shape and a 
good order-of-magnitude agreement on both axes. The 
agreement is surprisingly good taking into account that 
the aperture sizes are on the order of or larger, and 
thus our pinhole approximation is not justified. 

The more recent experiment of Ref.|5[was done with an 
array of apertures. In contrast to Ref. |j, this experiment 
should be well in the region where the linear approxima- 
tion is valid. However, a clearly nonlinear I— P curve 
was measured. The experimental results are also different 
for the two possible states of the weak link, the "H" and 
the "L" states, which have previously been identified as 
two nearly degenerate textural states with different R^'^^'s 
0. Yet, our numerical calculations confirm that loiP) 
curves in Fig.[51remain practically unchanged for all R 
This can also be seen in the adiabatic model as follows. 
The spin splitting eiSn-{4') ~ ii5n+{4> — V') has an es- 
sential effect on the critical current so that if can even 
vanish In ^Oj however, the relative phase shift is 

unimportant since the phase runs through all values. 

We conclude that the nonlinearity and the texture- 
dependence in Ref. are either large-aperture effects, 
and/or result form something other than MAR. In the 
following we demonstrate that at least part of the H-L 
difference can be accounted for by the anisotextural ef- 
fect, which exists in array-type weak links of ■^He-B. 

The anisotextural model. — As discussed above, the 
current-phase relation /(</>) depends on the rotation ma- 



3 



trices R ''^ . The basic idea of the anisotextural effect is 
that R''"^ are not fixed, but tend to move toward their 
(/>-dependent equilibrium configuration. This effect was 
previously used to explain the so-called tt states, where 
a local minimum of energy appears at (/> tt Below 
we generalize this theory to the dynamical case where 

(f){t) — LOjt. 

A simple but still realistic model for the anisotextural 
effect in an array of apertures is based on the energy 
functional Q 



(5) 



Here Fj is the Josephson coupling energy and the sec- 
ond term is due to the bending of the rotation axis n 
which parametrizes the matrix R(n, The rotation 
angle is fixed to 9l ^ 104° by the bulk dipole-dipole in- 
teraction Due to the geometry of the experiment 
0, we assume n to be fixed parallel to the wall nor- 
mal z on one side of the array. On the other side n 
makes an angle ri{r) with z, which depends on the dis- 
tance r from the center of the array. We choose a cutoff 
at r = To, below which 77 (r) is assumed constant (see inset 
in Fig. 12)). The Fj term depends on the angle 770 = 'ri{ro) 
near the junction, and the bending term is nonzero if 
770 differs from 7700 = preferred by walls and other 
textural interactions. At temperatures near Tc the en- 
ergy Fj can be approximated by Fj = —Ej cos (f> where 
Ej = aR'^,R'^,+ mUKa. + R\yRly)- Wc choosc to 



fix the r side (fi'' || z), and on the 1 side it is convenient 
to take < 77tx3 < 7r/2. The H and L states are then 
identified as the states with n'' = ±z, respectively 
Linearizing around 77 « 7700 , we find Ej « E^ — Jsp'&i'i'o)^ 
where 7? = 77 - 7700, and Jspirjoo) = [±(20a + 5/3) sin 27700 
-1-30/3 sin 7700] /16. The coefficients a, (3 and K are calcu- 
lated in Ref. |3 

The dynamics is determined by the Leggett equations 
S = -5F/Se and 9 = HollS/x Ell- Here the spin S 
and spin-rotation angle 9 = are conjugate vari- 

ables. 7g is the gyromagnetic ratio and x the spin sus- 
ceptibility. [The Leggett equations are canonical equa- 
tions for a Hamiltonian where the spin energy density 
fJ-olgS"^ f^X is added to F in Eq. ©J. The dynamical 
equations reduce to a wave equation in the bulk and Fj 
determines a boundary condition at the junction. The 
solution at constant P consists of waves that propagate 
radially out of the junction. This implies that the super- 
current Is = {2m^/K)d^Fj{rjQ, (j)) has a nonzero average, 
the anisotextural current 
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h 47 H- {ojjrQ/cf 



(6) 



Here 



y^2/jfl7ji^y(5x) is the spin-wave velocity and 

7 = birKro. We estimate tq = y/ A/tt, where A is the 
surface area of the array. We also take into account that 
the actual geometry around the junction differs from a 
half-space and use b = 0.31 instead of unity j3|. 
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FIG. 3: Comparison of experimental j3 and theoreti- 
cal (solid lines) current-pressure relations, (a) The L state 
current is compared with theoretical 7o (0Jl with a single ad- 
justable parameter a — 1.6 coming from the collision rate F. 
(Alternatively, one can neglect the gap suppression and use 
a = 3.2.) (b) The difference between the H and L state cur- 
rents is compared with the difference of Jai which is arbi- 
trarily multiplied by factor 3. The inset depicts the Josephson 
array, cut-off radius ro, and radiation of spin waves. 



The crucial features of the anisotextural current /ai 
© are the dependence on the textural configuration via 
JspiVca), and the nonlinear dependence on pressure via 
LUj. The origin of /ai is the oscillating part in Eq. 
which was here approximated by the if term alone. The 
total dc current I^c is the sum of /ai and the constant 
component due to MAR, /dc = lo + /ai • 

Comparison to experiments. — It turns out that the L 
state current, which is rather linear can almost com- 
pletely be attributed to MAR. The comparison of the L 
state data with Iq is presented in Fig.|2Ia). There is only 
one adjustable parameter a ~ 1, which is associated with 
the uncertainty of the relaxation rate F. A likely rea- 
son for the differences is that the theory uses the pinhole 
approximation, which is not strictly valid for the exper- 
imental apertures (100 nmxlOO nm squares in a wall of 
thickness 50 nm are not small compared to ^0 = 77 nm). 

Since both the H and L state data are assumed to con- 
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tain the same contribution Iq , it is convenient to subtract 
the H and L state data taken at equal temperatures. In 
this way one may directly compare the texture-dependent 
parts with /ai , regardless of any uncertainty that may be 
present in Iq. This is done in Fig. |3{b). The theoreti- 
cal result corresponds to the difference between Eq. © 
calculated for the two states. There are no adjustable pa- 
rameters, since the only free parameter is the textural an- 
gle rjoo , which was previously found to be approximately 
O.Stt based on static properties Q- The anisotextural 
model, as presented above, can explain roughly one third 
of the observed H-L difference. [The theoretical current 
in Fig. Olb) is arbitrarily multiplied by factor 3 in order 
to make it better visible.] Also, the curvature in the the- 
oretical lines is at higher biases than in the experimental 
results. There are several possible sources for the dif- 
ferences, in particular the oversimplifications used in our 
anisotextural model. Unfortunately, it would be very de- 
manding to improve upon the pinhole approximation, or 
to calculate the texture and propagation of spin waves in 
the complicated geometry of the experiment There 
is also uncertainty in the experimental parameters, for 
example in the diameter of the apertures, which appears 
in its fourth power in /ai. As a result, the true reason for 
the factor-of-three discrepancy remains open. 

Conclusion. — We have presented a theory of dissipa- 
tive currents in '^He weak links. It shows that the non- 
linearities in the measurements of Refs. and have 
different origins, and both can semi-quantitatively be ex- 
plained by natural extensions of existing theories. The 
extension of static anisotextural phenomena to dynam- 
ics gives further support for the theory, and provides the 
energy loss mechanism that is required for a weak link 
to become trapped in a tt state, as seen in experiments 
Q. The anisotextural phenomena are sensitive to the ex- 
perimental cell and magnetic field, for example, and thus 
can be tested in detail in future experiments. Also the 
oscillating components and in Fig. ^as well as the 
spin waves might be observable in experiments. 

Appendix: Calculation of currents. — Consider a pin- 
hole with open area Aq. Assume the 1-side chemical po- 
tential to be shifted by U with respect to the r side, 
and take the z axis to point from 1 to r. We define 



Gn in Eq. Q as G„ = ^m3VFN{0)Ao, where 7V(0) is 
the single-spin density of states in the normal state and 
vp the Fermi velocity The mass current may then 
be written as 1(1) = Ga{kJ{k,t)) f,^^^ where (■■•)fe^>o 
~ /fc >o(*^^k/^^) ' ' ' denotes an average over the Fermi- 
surface points k. Since I{t) is periodic with period 
Tj = we expand /(k,<) = E,T=-oo ^n(k)e'""^* 

such that /„(k) = /*„(k). Using the 7 matrices of 
Ref. to expand the Keldysh function, the amplitudes 
for n ^ may be written as /„(k) — TrC{2[/(5„o — 

E™=o^/d£E"'"^"(k,e,c/) - z:r"^'"(-k,£,-^)]} 

where F[f (k,e,C/) = - 7'^*(e)i^> - 

t/)7^'(e)][£[,(e,C/)]t andP™(e,[/) = n^Jo' 7^^(e+(2p+ 
2)f/)7^^(e + {2p + 1)U) where ij = l,r, C = 1, and k 
dependences are omitted for clarity. (For spin current, 
replace TO3 h/2 and Since we are considering 

a pinhole, the 1 and r order parameters may be assumed 
to be unaffected by the weak link and to be in equilib- 
rium. Thus the distribution functions are of the form 
x^(e) = h'{e) [1 -7^*(e) l^'^^)] ^'^^) = ^ " ^/(e). 
The conjugation symbol " ~ " is defined as in Ref. jf^ 
and the amplitudes satisfy 7^ = [7^]^- 

Above the 7-'^''^'s refer to the coherence functions in- 
side the junction. They are obtained by integrating the 
Riccati equation i/ivp • V7^ = -2e^7^ "I^AoIo ~ ^0 
on several trajectories toward the wall on both sides, 
where = e + i(^r/2), Aq = Ao • ctIctj, Ao(k, z) = 
(Aiifca;, Ayfcy, A^fc^), and A||(2:), A_l(2;) are calculated 
with the "ROM" model of a diffusive surface (cf. Ref. 
3]). For a symmetrical junction, the 1 and r solutions at 
the junction satisfy 7^'(-k) = -[7^'(k)]'^. The differ- 
ent spin rotations C/'"'^ = cxp(— i^^^n''"^ • g_/2) arc taken 
into account by writing 7^* = II^1q^\U^V ^ for i = 1, r. 

The current amplitudes of Eq. iQJ in Fig.[5]are obtained 
with /o = G'„(^^Jo(k))j^>o, II = 2GnIm(A:^/„(k))^^>o, 
I,'^ = 2Gn Re(fcz/„(k))j, ^p. These results have been de- 
rived for time- independent R''"^, but at least for U <C HTq 
we may expect them to be valid even if the textures vary 
on time scale Tj. 
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